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Abstract 

Electroweak baryogenesis depends on the profile of the bubble wall created in 
the first-order phase transition. It is pointed out that CP violation in the Higgs 
sector of the MSSM could become large enough to explain the baryon asymmetry. 
We confirm this by solving the equations of motion for the Higgs fields with the 
effective potential at the transition temperature. That is, we present an example 
such that the transitional CP violation is realized and show the possibility that the 
baryon asymmetry of the universe may be produced, if marginally, by the r lepton 
interacting with the wall, when an explicit CP breaking in the Higgs sector, which 
is consistent with experimental bounds, is induced at the phase transition. 
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1 Introduction 



BaryogenesisfQ at the electroweak phase transition (EWPT) is an attractive mechanism 
to explain the baryon asymmetry of the universe (BAU) because the mechanism depends 
on parameters that will be tested sooner or later on the earth. For the mechanism to be 
viable, however, some extension of the standard model is required to guarantee first-order 
EWPT and to incorporate other sources of CP breaking than the CKM phase. 

Among possible extensions of the standard model, the minimal supersymmetric stan- 
dard model (MSSM) may be a well-motivated one. The MSSM could cause a strongly 
first-order phase transition when one of the stops is light [0, [3| and has some sources of CP 
breaking, though they are constrained by measurements of the neutron electric dipole mo- 
ment (EDM) j|, |5|. On the other hand, another source of CP violation is relative phase 
6 of the vacuum expectation values (VEVs) of the two Higgs doublets[|]. The Higgs 
VEVs including the phase, which characterize an expanding bubble wall created at the 
first-order EWPT, vary spatially. Such CP violation affects propagation of quarks and 
leptons through their Yukawa couplings and that of charginos, neutralinos and sfermions 
through their mass matrices, so that weak hypercharge is carried by these particles into 
the symmetric phase region, where the hypercharge will be turned into baryon number 
by sphaleron processes. 

In previous papers^ §, we attempted to determine the profile of the bubble wall by 
solving equations of motion for an effective potential V e g at the transition temperature Tc 
in the two-Higgs-doublet model under some discrete symmetry. For a phenomenological 
set of the effective parameters, we presented a solution such that the CP-violating phase 
9 spontaneously generated becomes as large as 0(1) around the wall, while it completely 
vanishes in the broken and symmetric vacua. We shall refer to this mechanism as tran- 
sitional CP violation. This solution yields the hypercharge flux by the quark or lepton 
transport to generate the BAU. We also showed f| that a possible explicit CP break- 
ing in the Higgs sector at Tc does nonperturbatively resolve the degeneracy between the 
CP-conjugate pair of the bubbles and leave a certain amount of BAU after the EWPT. 
Although several necessary conditions for the transitional CP violation were found in the 
MSSM[|l(|][], we must solve the equations for the wall profile in order to confirm that such 
a mechanism works in practice. 

In the present article, we investigate possibility of the transitional CP violation in 
the MSSM. First of all in §2, we study the effective potential at Tc following a method 
proposed by one of the authors(K. F.) ||12|| . It is approximated by polynomials in the Higgs 
fields whose coefficients are given by the effective parameters (1^, [13], [14]]. In particular, 
we also point out that an explicit CP breaking in the Higgs sector, which is induced 
through loop corrections from the SUSY particles, could be enhanced at the EWPT. 
Employing the approximated effective potential, the equations of motion for classical 
field configurations are derived to find a CP-violating profile of a bubble wall in §3. In §4, 
we give one of numerical examples of the effective parameters which, at zero temperature, 
produce mass spectrum consistent with present observations. The global structure of 
the effective potential suggests that two degenerate minima of it are connected by a 
path with almost constant tan/5, where tan/? is the ratio of the expectation values of 

1 As for studies including nonperturbative effects by use of 3d effective theory, sec jy] 
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the two Higgs doublets^. This justifies to neglect spatial dependence of tan (3 so that 
the number of unknown functions is reduced. Though the bubble wall is rather broad, 
transitionally CP- violating solutions by the effective parameters may produce the BAU 
by the r lepton transport, if marginally, when the explicit CP breaking in the Higgs sector 
is induced from complex parameters consistent with the neutron EDM measurement. A 
summary and discussions are given in §5. Formulae for the effective parameters are given 
in Appendix A. 

2 Effective Potential and Transitional CP Violation 



Let us parameterize the Higgs doublets of the MSSM as 
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where 9 = 9\ — 02- The most general potential for them at the tree level, which is gauge 
invariant and renormalizable, is given as 
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Here the coefficients are as follows: 
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where g 2 (i) is the SU(2)(U(1)) gauge coupling and fi is the coefficient of the Higgs 
quadratic interaction in the superpotential. The mass squared parameters m 2 , m\ and \xB 
come from the soft SUSY-breaking terms so that they are arbitrary at this level. m\ could 
be complex but its phase can be eliminated by rephasing the fields when A5 = Xq = Xj = 0. 
We adopt the convention in which this m\ is real and positive. 



2.1 Effective potential 

The effective potential up to the one-loop terms of radiative and finite-temperature cor- 
rections is 

V cS = V + V 1 (p l ,0) + V 1 (p i ,T), (2.5) 

where 
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2 It is also reported that the spatial dependence of (3 is very small: A/3 ~ O(10 3 )[[L5| 
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Here we consider the contributions of gauge bosons, top quark, top squarks (t), charginos(x ± ) 
and neutralinos The n, counts the degrees of freedom of each species including its 

statistics, that is, rij > 0(rij < 0) for bosons (fermions). The rrij, which is a function 
of the Higgs background (pi,9), is the mass of the particle. The mass matrices of the 
charginos, neutralinos and stops are given in Appendix^. 

We assume Mi = M 2 for the gaugino mass parameters in the mass matrices, which 
simplifies the neutralino contributions so as to be proportional to the chargino contribu- 
tions. The coefficients in Vq are replaced by the corresponding effective parameters in 
Veg. Such effective parameters that are relevant to CP violation in the Higgs sector are 
defined by 
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where A x (Aj) implies the one-loop corrections from charginos and neutralinos (stops) 
given in Appendix. Except for the light stop (t\) contributions, Taylor expansion of the 
finite-temperature integrals around pi = is valid. As for the light stop, we employ the 



high-temperature expansion ||16|| to obtain 
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where a|(p) = m|(p)/T 2 and A_ = 0.1764974. 

Finally, the effective potential up to the p 4 -terms is expressed as 
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3 We have used the DR-scheme to renormalize T4ff. 
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where rhi and \ are the effective parameters, some of which are given above. The 
^-dependent p 3 -terms come from the light stop contributions by the high-temperature 
expansion [[10[|, while the weak gauge bosons contribute to ^-independent p 3 -terms. Note 
also that non-zero values of A5, \ and A 7 are induced. 

To calculate the effective potential, we give v = 246GeV and tan (3 at zero temper- 
ature as an input. Once m\ at the tree level is given, the other mass parameters in the 
tree-level potential is fixed by minimizing the effective potential at T = 0||12|. At the 
same time, the mass of the lighter (heavier) CP-even Higgs scalar (m#) and that of 
the pseudoscalar uia are evaluated from the second derivatives of the effective potential 
at the minimum. When the EWPT is of first order, the transition temperature Tq is de- 
termined as temperature at which (pi,p 2 ) = (0,0) and (pi,p 2 ) = (v cos /3, v sin (5) become 
degenerate minima of the effective potential, where the symmetry-breaking minimum is 



searched for by a numerical method explained in [FJ|. This can easily be extended to 
the case with CP violation, in which 9 acquires nonzero value and the Higgs scalars and 
pseudoscalar mix to produce new mass eigenstates[|T^, |H| . 

2.2 Spontaneous CP violation 

For simplicity, we assume that all the parameters of V e s are real, that is, no explicit CP 
breaking. Then at zero temperature CP symmetry can be spontaneously violated but 
the lightest scalar becomes too light to satisfy experimental lower bound |1S||. 

Now we argue that if the EWPT is of first order, CP can be violated for broader range 
of acceptable parameters than at T = 0. Note that 

V e ff{pi,0i) = F(pi,p 2 ) [cos 8 — G(pi,p2)} 2 + ^-independent terms, (2-15) 

where 

F(pi,P2) = yP?p^-2(S 2 p?p 2 + C 2 p 1 pl), (2.16) 

n( \ 2m| + \p\ + \ lP \ + 2(B lPl + C lP2 ) 

G(Pl ' P2) = 2X 5PlP2 -8(B 2 p 1 + C 2 p 2 ) • (2 ' 17) 

Recall that spontaneous CP violation occurs at T = {B{ — Cj = 0), only if 

F{p u p 2 ) > and - 1 < G(pi, p 2 ) < 1, (2.18) 

for (pi, p 2 ) = (v cos /3o, vq sin p ). These conditions strictly restrict the range of m| so that 
a very light scalar inevitably results, since fh\ must be the same order as \§p\ + A7P2. At 
the first-order EWPT, (pi, p 2 ) varies from {y cos f3, v sin (3) to (0, 0) between the broken and 
symmetric phase regions. Then the conditions ( |2.18| ) could be satisfied for broader range 



of parameters, which are acceptable in contrast to the zero-temperature case, since the 
effective parameters receive finite-temperature corrections and (pi, p 2 ) is not specified to a 
single point. In this case, there exists a CP-violating local minimum in the transient region 
and the bubble wall profile, which is represented by the classical Higgs configuration, could 
have nontrivial 9 near such a minimum. 
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Even when P(pi,p 2 ) < for some (pi,p 2 ), the bubble wall profile could acquire 
nontrivial CP phase as long as —1 < G(pi, p 2 ) < 1 is fulfilled at that point. In this 
case, V e f{ for the fixed (pi, p 2 ) is convex upwards and its peak is located at 9 G (0, n). If 
G(pi, P2) varies between positive and negative values along the wall profile, the minimum 
of V e fi for a fixed (pi, p 2 ) travels between 9 = n and 9 = as pi varies from the symmetric 
phase to the broken phase. Then there exists a CP- violating saddle point of V e g in the 
transient region. 

In any case, to realize —1 < G(pi,p 2 ) < 1 at some (p 1; p 2 ), \fh\\ at T c must become 
much smaller than tree-level m\ by large negative corrections. Since the correction from 
the charginos (stops) to m| is proportional to pM 2 (pA t ) as shown in Appendix, pM 2 < 
and/or pA t < are required for transitional CP violation to occur. Once these parameters 
are provided, m| which leads to — 1 < G(pi, p 2 ) < 1 is restricted in a rather narrow range. 
Then the mass of the Higgs scalars and pseudoscalar are almost uniquely determined. 
Whether a transitionally CP- violating bubble wall is realized or not is determined finally 
by solving the equations of motion with the effective potential at Tc- 



2.3 Explicit CP breaking 

In the discussions before, all the parameters of V e s are assumed to be real. However, an 
explicit CP breaking is necessary to avoid the complete cancellation in the net baryon 
number due to the symmetry 9 <-> —9 between the wall profiles. 

The origins of the explicit CP breaking in V e s are those in the complex p-parameter 
and SUSY-breaking parameters which depend on the phases, a\ = arg(pMi) = arg(pM 2 ) 
contributed from charginos and neutralinos and a 2 = axg(pA* t ) from stops, as given in 
Appendix. From these we have 
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where A^ ^ denotes the corrections in the case when all the tree-level parameters are real. 

The magnitudes of a± and a 2 will be bounded by experiments to be 0(< 1CT 3 ) or 
so[|J. These CP-breaking phases can be gathered on m\ by rephasing the Higgs fields 
when the contributions of the stops are very small compared with those of charginos and 
neutralinos. After rephasing, m| is given in the form of 
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so that the CP-breaking phase S is given by 
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•C mj, which is often 



This suggests a very interesting possibility that, if mj + A^'mi 
the case for the transitional CP violation, then a somewhat large \S\ is induced even if 
sinai ~ 0. 
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Hereafter, we put the explicit CP breaking at Tc by replacing m\cos6 in V e s (|2.14[ ) 
by (l/2)(e l( - 5+0 ^ml + h.c), but consider the case of small 5 ~ O(10 -3 ). 

3 CP- Violating Solutions 

3.1 Equations of motion 

We are now interested in classical solutions of the bubble wall. If the phase transition 
proceeds calmly, it will be valid to expect that the bubble wall grows keeping the profile 
of the critical bubble, which is determined by static equations of motion. Further, when 
the bubble is spherically symmetric and is sufficiently macroscopic so that it is regarded 
as a planar object, the system may be reduced to an effectively one-dimensional one. 
Regarding the bubble wall as a static planar object, the equations of motion are given 
byi 



d 2 pi(z) 
dz 2 




Pi( z ) ( "J— ) " — = °> 

= 0, (3.1) 

where z is the coordinate perpendicular to the wall. Furthermore, gauge configurations of 
the pure-gauge type with no dynamical freedom in 1+1- dimensions are expected to give 
the lowest energy of the system. Then, it is convenient to fix the gauge in such a way 
that the gauge fields disappear by imposing the condition 

pl{z)^- L +pl{z)^ 1 = Q. (3.2) 
dz dz 

The total energy density of the bubble wall per unit area is given by 



£ 



3.2 Kink ansatz 




dfiV 2 ( dOiV 
dz I ^\dz) 



+ Veff(pi,P2,0)>. (3.3) 



In accord with the first-order phase transition, we adopt an ansatz that the kink config- 
uration connecting the vacua is a solution to the equations of motion. Using a dimen- 
sionless variable y = (1/2) (1 — tanh(a2;)), we put p(y) = (pi/v)/ cos/? = (p 2 /v)/ sin/5 and 
9(y) = 0\/ sin 2 /? = — 9 2 / cos 2 (3 under the gauge-fixing condition (|3~2|). Here a is taken to 
be the inverse wall thickness. The equations of motion to give p(y) and 8(y) are read as 

4,(1 -v)%(l -,)*)-4cos^sm 2 Ml-,)V(? 2 - ° W(P 



dy^ dy' ^dy' dp 

4cos /3sm py{l-y) — [y{l-y)p — ) = — , (3.4) 
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where W(p(y),9(y)) = V c s(pi,9)/a 2 v 2 is dimensionless. 

The parameterization of W convenient to the kink ansatz in the presence of 5 is as 
follows: 
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Here ~ 0(5) is the boundary value of 9 in the broken vacuum given later. The 
parameters in W(p, 9) and those in V e s(pi, 9) are related as follows: 
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together with relations from the kink ansatz: 

m\ cos 2 13 + m 2 . sin 2 (3 = 2a 2 (2 + 6cos(6» 6 + 5)), 

Ai cos /3 + A 2 sin 4 /3 + 2(A 3 + A 4 ) cos 2 /3 sin 2 /3 

= {a/v) 2 8{2 + ccos9 b + {d/A)cos29 b ), 

A cos 3 /3 + (B cos(3 + C sin/3) cos /3 sin /3 + Dsin 3 /3 

= (a 2 /v)(4 + ecos6» f , + /cos26' 6 ). 



(3.7) 



Once the parameter set (b, c, e, d, f) is given, the stability of the two vacua and the con- 
dition W(p, 9) > in the region of < p < oo and < 9 < 2tt have to be checked. As it 
should be, in the case of 5 = 0, the kink solution (p — 1 — y = p^nfe, = 0) satisfies the 
equations of motion. 

That dW/d9\ p= i = gives the boundary value 9 b determined from 



tan 6*6 = 



b cos 



b sin 5 

5 + c - e + (d - 4/) cos 9 b ' 



(3.8) 



and dW /d9\ p ~o = does the boundary value of the symmetric vacuum as 9 S = —5, ir — 5. 
Needless to say, p b — 1 and p s = 0. At first sight, any 9 S might be allowed because of 
p s = 0. But the energy density of the bubble wall diverges for the other 9 S . 

In the case of 5 = 0, suppose that we obtain a nontrivial solution (p ^ Pkink, MO) 
with the maximum of 9 > 0. Let us denote this solution as (p°, 9°). As explained before, 
we have the CP-conjugate partner of it, (p°, —9°). Without loss of generality we choose 
5 > because of the symmetry (5,9) < — >■ (—5, —9). For 5 > small enough, we would 
have in general three types of solutions, the positive-6 1 solution (p + ,9 + ), the negative-6* 
solution (p~,Q~) and the small-0 solution (p s ,9 s ) such that, as 5 — > 0, the positive-6 1 
solution tends to (p°,9°), the negative-6 1 solution to (p°, —9°) and the small-6 1 solution to 
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Table 1: Parameters adopted for the numerical calculation. 



^0 


tan p Q 


m\ p 


A t 


M 2 = Mi mi L m- tR 


246 GeV 


6 


8110 GeV" -500 GeV 


60 GeV 


500 GeV 400 GeV 



the kink solution (pkink, = 0). The small-6 1 solution has no CP-conjugate partner. For a 
larger 5, either of the positive-6 1 solution or the negative-6 1 one, which has the lower energy 
for a small 5, survives 0. 

The energy density of the bubble wall is now expressed as follows, where that of the 
kink solution, Ekink = av 2 /3, is chosen as the standard: 
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dy, 



(3.9) 



Let us denote as AS + , AS and AS S respectively for the positive-^, negative-^ and small-^ 
solutions. 



4 A Numerical Example of Transitional CP Violation 

4.1 Numerical results of MSSM calculation 

We give one of numerical examples at the EWPT and the effective parameters in the 
MSSM. The parameters adopted are listed in Table [[]. As noted in §2.2, pM 2 < and 
pA t < are favored for transitional CP violation. Once these parameters are given, the 
tree- level m 2 is tuned to satisfy the conditions for transitional CP violation. Since the 
magnitudes of the effective couplings at Tq ~ 100 GeV are always 10~ 3 — 10~ 2 when 
the mass parameters are taken to be between weak scale and a few TeV, an appropriate 
value of nil is about O(10 4 ) GeV 2 . Although it is reported that a small value of tan/9 is 



favored for a strong EWPT satisfying v/Tc > 1 with an acceptable Higgs mass|L9]], we take 
tan Po = 6. One of reasons why we take tan Po = 6 is that for the other parameters listed in 
Table [|, a larger tan p yields a heavier Higgs scalar when one includes contributions from 



charginos and neutralinos, which were often ignored in literatures |T^|. Another reason is 
that for a large tan/3o, it is rather easy to find a stable nontrivial wall profile, since the 
second kinetic term in ( p.9|) ~ cos 2 /? while W(p, 6) ~ cos/5 as in (|3.6|) , so that AS < is 
realized. 

The particle spectrum from this set of the parameters are given in Table |2|. The 
transition temperature Tq and the VEVs at Tc are given by 

T c = 93.4 GeV, v = 129.17 GeV, tan p = 7.292, (4.1) 

so that v/Tq — 1.38, which guarantees that the BAU created at the EWPT is not washed 
out after the transition. The values of the effective parameters near Tc are presented in 



8 



Table 2: Masses of the lighter Higgs scalar (rrih), the Higgs pseudoscalar (m^), the heavier 
Higgs scalar (m#), the lighter top squark (mjj, the lighter chargino (m^) and the lightest 
neutralino (to % o). 

rrih rnj itlh m xt m *i 

82.28 GeV 117.9 GeV 124.0 GeV 167.8 GeV 457.6 GeV 449.8 GeV 



Table 3: The effective parameters near the transition temperature Tc = 93.4 GeV. 



T(GeV) 


(ml) cS (GeV') 




A 5 




A 6 




A 7 




93.3 


-402.375 


4.95 


3562 x 10" 


-3 


-1.09172 x 10" 


-2 


7.65974 x 10- 


-3 


93.4 


-401.720 


4.95 


3684 x 10" 


-3 


-1.09197 x 10- 


-2 


7.65912 x 10- 


-3 


93.5 


-401.065 


4.95 


3806 x 10- 


-3 


-1.09222 x 10- 


-2 


7.65850 x 10" 


-3 


93.6 


-400.410 


4.95 


3928 x 10- 


-3 


-1.09247 x 10" 


-2 


7.65787 x 10- 


-3 



Table |3| and 

B C 

— = 1.5741 x 10~ 3 , — = 3.2539 x 10~ 2 . (4.2) 

Tc T c 

The global structure of V e s at Tc is depicted in Fig. [I]. The T dependence of the effective 
parameters are plotted in Figs. and |3|. 

By closely investigating the saddle point between the two degenerate minima, we find 
that the height of the saddle point is V max = 3.6516 x 10 9 GeV 4 , from which the inverse 
wall thickness a is given by 



a = v/8Kiax = 13.23 GeV. (4.3) 

v 

4.2 Transitional CP violation and BAU 

From the numerical results above, we fix the parameters at 7^=93.4 GeV as -y=129.2 
GeV, a=13.23 GeV, tan/3=7.29 and (b = -0.308, c = 0.048, d = -0.009, e = -0.3,/ = 
—0.0003). For 5=0.001 and 0.002 respectively, a pair of solutions (p + (y) , 6 + (y)) and 
(p~(y), 9~(y)) with 9 S = ir — 5 are obtained, which are shown in Fig. f| of the W1-W2 plain 
on the contour plot of W(wi,W2), where W\ = pcos(6 — 6b) and W2 = psin(# — Bb) 0- 
Bubble formation rate The bubble formation probability of the each solution is 

4 The solution for 0" = ^ with V e g = \§ 2 {<\> - vf is given by <f> = |(1 + tanhaz) with a 2 = |w 2 . 

This leads to V max = V e s{<j> = f ) = ±v 4 = ^f. 

5 For b < 0(> 0) it is easy to show that the exit gate toward w\ < 0(lcft) from the basin of the 
symmetric vacuum at = 0, W2 = 0) is wider (narrower) than that toward w\ > (right) in Fig. |^. So, 
b < 0(> 0) favors 9 S = n — S(= —5). Usually we have found no convergent solutions with 9 S = —5(tt — S) 
for b < 0(> 0) in the relaxation method. For b < with tan/3 < 3, no solutions have been obtained. For 
the smaller values of tan/3 there is a wide parameter range to admit solutions for b > 0. However, the 
parameters of them, in particular (c, d, e), are not compatible with the MSSM calculation. 
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Figure 1: Contour plot of V c g(pi;Tc) in the unit of GeV 4 , where p\ 
v sin (3 in the unit of GeV. 
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Figure 2: Behaviors of the effective parameters (fhf) eS and As. The dashed lines represent 
the contributions from the charginos and neutralinos, the dotted lines represent those from 
the stops and the solid lines are the sum of the tree-level and one-loop contributions. 
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Figure 3: Behaviors of the effective parameters Xq and A7. The dashed lines represent the 
contributions from the charginos and neutralinos, the dotted lines represent those from 
the stops and the solid lines are the sum of the tree-level and one-loop contributions. 
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Figure 4: Contour plot of W{wx,w^) together with the positive-^ solution for S = 0.001 
(broken curve). 



proportional to 

N ± = exp ( |g J , (4.4) 



where R cr it is the radius of the critical bubble at T c and is given by ^j3F crit /4:irav 2 and 
the free energy of the bubble F crit is estimated as 145Tc. The explicit CP breaking 5^0, 
even when very small, nonperturbatively breaks the degeneracy of the CP-conjugate pair 
of the bubbles as mentioned before. The formation rate is 

N-/N + = 0.601 for 5 = 0.001, 

= 0.361 for 5 = 0.002. (4.5) 

Chiral charge flux Once a set of solutions are given, we calculate the reflection 
coefficient AR(£,pl/o) of the each solution, where pi is the incident momentum of the 
particle perpendicular to the wall and £ = mc/fl, m c being the mass at T c . That is, 
mc = m fc ( r )t>cos/?/(voCOS/9o) for the bottom quark (r lepton). Important here is that the 
bottom quark and the r lepton interact with the wall magnitude p(y) with the coupling 
sin 2 /3 while the top quark does with cos 2 /?. Therefore the top quark passes the wall 
almost freely for such a large tan j3 and is completely negligible for the BAU. 
Then the chiral charge flux through the wall is given by|[20|| 

r Ql ~ Qr 



47T 2 7 



f°° dp L r d p T p T [f s -f b }AR(^) } (4.6) 
Jm c Jo 1 J a 



where 



fs = P± 1 fb = PL 1 U 7) 

E e^E-u PL )/T c + i ' ■> E el{E+u ^prz^r )/Tc + 1 > v • > 

are the statistical factors in the symmetric and broken phases respectively, P t is the 



incident momentum of the particle parallel to the wall, E = Jp L + P j., u is the wall 



velocity and 7 = l/yl-M 2 . 
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Figure 5: Contour plot of F$ et /((Q L - Q R )T%u 



The net contribution to BAU is [] 



rpnet 

Q 



iV+F+ - N-F Q 
N+ + N- 



(4- 



Fig. | shows FQ et /((Q L - Q R )T^u) for 5 = 0.001 and u = 0.1 as a contour plot on 
the a/Tc — mc/Tc plain. The point B + (B_) is the b quark contribution if m=4. 4(4.1) 
GeV, and the point T is the r lepton contribution. Fq et /((Q L — Q R )Tjjii) amounts to 
4 x 10~ 7 ,3.4 x 10" 7 ,4 x 10~ 8 at B+,B_, T respectively. The top of the contour hill is 
6.3 x 10~ 3 at log 10 (a/Tc) = 0.88 and log 10 {m c /T c ) = 0.29. As is well known, a thin wall 
produces a large chiral charge flux. For 5 = 0.002, F"q 1 increases as a whole by a factor of 
about 1.7. As u increases to 0.9, the contour at the point £>_, say, decreases by a factor 
of about 0.4. 

From Fp et we can estimate the BAU in the transport scenario as 



^ ~ io- 7 



rpnet 
X — — X 

u 



m2 i 
1 C 



(4.9) 



where the entropy density is given by s = 2tt 2 ^T^/45 with ~ 100 and r is the 
transverse time within which the scattered particles are captured by the wall. If r is 
given by D/u with the diffusion constant D, it is estimated from D(b quark) ~ 1/Tc 
and D(t lepton) ~ (10 2 ~ 3 )/T C . Then 



Pb/s 
Pb/s 



< 



10 " (for b quark), 
10 -(io-i2) ( for T lepton ) 



(4.10) 



6 For the parameter set with b < 0, we have found no small-# solutions because the wider exit gate 
from the symmetric vacuum has the direction just opposite toward the broken one at (wi = 1,W2 = 0). 
Even in the case when the small solution exists, A£ s is negligibly small. 
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for 5 = 0.001 and u = 0.1. Thus our example may give a possibility to produce the BAU 
by the transitional CP violation of the r lepton, if marginally, in the MSSM, provided 
that 5 ~ O(10~ 3 ) is induced. 



5 Summary and Discussions 

We have shown the possibility of the transitional CP violation in the MSSM. For this 
scenario to be realized, some of the parameters in the theory are constrained: \lM-i < 
and/or jiA t < are required, m| should be O(10 4 ) GeV 4 , and tan(3 > 5 would be 
necessary for our choice of the mass parameters to have the lightest Higgs scalar with an 
acceptable mass and to have nontrivially CP-violating wall profile. These requirements 
inevitably relate the mass of the Higgs scalar to that of the pseudoscalar. 

Although the explicit CP breaking in the Higgs sector is severely limited from the 
neutron EDM at T = 0, that induced at Tq can be somewhat large. If the latter is 
5 ~ O(10 -3 ) and breaks the degeneracy between the pair of CP-conjugate bubbles, our 
numerical example could be able to produce the BAU~ O(10 -10 ) by the r lepton trans- 
port. The ^-dependent p 3 -terms essential for the transitional CP violation is induced only 
if one of the soft SUSY-breaking masses, mi R of the stop, almost vanishes. This implies 
that the mass of the lighter stop must satisfy < m t . 



Our result of a rather thick walls, 1/a ~ 10/Tc together with a large A9 = \9b — 9 s 



~ 7T 



may favor the diffusion scenario |2T|, although estimating various uncertain factors of this 
scenario is outside the scope of this article. 

In §2.2, we have pointed out another possibility of transitional CP violation associated 
with P(px,p 2 ) < 0. A wall profile of this type, if exists, will generate sizable BAU 
just as the example presented here, since 9 varies from to ir. As long as the high- 
temperature expansion of the light stop contributions is valid, the P 2 -term in ( |2.16|) 
is negligible compared to the As-term. Indeed we found a well qualified solutions for 
(d = —0.009, / = —0.0003), which is obtained from the parameter set studied here only 
by changing the sign of d. It is, however, difficult to have negative A 5 and |G(pi,p 2 )| < 1 
for an intermediate (pi,p 2 ) with an acceptable set of parameters when mi R = 0. This 
can be seen as follows. As shown in [10], Apri 2 / \ixA t ) and AjAs/(pyl t ) 4 are functions of 



T and when mi R = 0. These are plotted for tan/5 = 6 at T = 95 GeV in Fig. [fj. 
Contour plots of the chargino-neutralino contributions A x rn| and A x As for tan/3o = 6 
and T = 95 GeV are shown in Fig. |7|. In order for A5 to be negative, rather small |/i.M 2 | is 
required. Since | A X A 5 | is at most O(10~ 4 ) in the region where A X A 5 is negative, \fiA t \ must 
be smaller than 10 4 GeV 2 to have negative A5. As seen from Fig. 0, A. x ml > 1500 GeV 2 
in the region where A x As < 0. As long as we take larger than the weak scale, 
AjmJ ~ —0.1 • fiA t , so that A^m| should be larger than —1000 GeV 2 . Hence to satisfy 
|C(pi,p 2 )| < 1, the tree-level m\ must be taken to be at most 2500 GeV 2 . This small 
m\ inevitably produces too light Higgs bosons, which are inconsistent with the present 
bound rrih > 67.5 GeV and > 67.5 GeV[^2|. We have examined cases with larger 
tan fio, which yield heavier Higgs bosons, but find no region in the parameter space that 
both A5 < and |C(pi,p 2 )| < 1 are satisfied with acceptable Higgs masses. Although this 
interesting scenario of CP violation is far from realizable in the MSSM, a general two- 
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Figure 6: Ap7%/(jiA t ) and A t ~A 5 /(Mt) 2 for tan ft = 6 at T = 95 GeV. 
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Figure 7: Contour plots of A x m| and A X X 5 at T = 95 GeV. 

Higgs-doublet model will allow a wall profile with A5 < 0, since it has nonzero tree-level 
A5 as a free parameter and we confirmed the existence of such a solution. 
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A Formulas for the Effective Parameters 

In this appendix, we summarize the contributions from the charginos, neutralinos and 
stops to the effective parameters. As for the formulas to evaluate the finite-temperature 
Feynman integrals, refer to Appendix of [[Uj . 
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The mass matrices of the charginos and neutralinos are given by 

/ M 2 -%(hfi ' 



W2 _ 

7^ 2 



// 
/i 



(A.l) 



(A.2) 



7i^ e 
-p 

( M 2 
Mi 
-f#2Pi 

respectively. Here M x and M 2 are the gaugino mass parameters. The contributions from 
the charginos are given by 

A x± m 2 = 2g 2 2f iM 2 -i [ A 1 (k)A^k), (A.3) 

J k 

A X± A 5 = -2g^M 2 ) 2 -i f A?(*)A£(fc), (A.4) 
A X± A 6 = A x± A 7 = 2^ 4 /iM 2 -2 [ k 2 Al(k)Al(k), (A.5) 

and those from the neutralinos are 

2% [ U^A^A^k) + gliiM.A^A^k)] , (A.6) 

-2i j \glnM 2 & x {h) + gf^A^k)] 2 A 2 (k), (A.7) 

Jk L J 

2* / A; 2 k 4 M 2 A 2 (A;) + <7 2 2 <? 2 (M 2 + M 1 )A 1 (A;)A 2 (A;) + <? 4 MiA 2 (A;)l fiA^k) 

Jk L 

A x oA 7 , 



A x om 2 
A x oA5 
A x o A@ 



where 



Ai(fc) 



A; 2 - iMol 2 ' 



A 2 (fc) = 



1 



A' 2 



I Mi 



12' 



1 



A; 2 — |/i| 2 



(A.8) 
(A.9) 



Note that M 2 , Mi and /i are complex. 

In the special case of M 2 = M l5 the corrections to the effective parameters are reduced 

to 



A x o m\ 



A x oA 5 = 



2 

COS 2 0w Jk 



COS 4 



(# 2 f^ / A?(*)Aj(fc), 

J k 



A x oAg 



A x o A7 



at 



-nM 2 -i [ k 2 Aj(k)Al(k). 

J k 



x " "cos 4 ^' 

In this case, the contributions of charginos relate to those of neutralinos as 



A x om 2 



1 



-A x ±m 2 3 , 



cos 2 9 W 
A x oA 5 = — — r-j A X ±A 5 , 



cos 4 6w 
A x oAg = A x oAy = 



A Y ±A 6 = — — — A Y ±A 7 . 
x cos 4 9 x 



cos 4 6 W 
15 



w 



(A.10) 
(All) 
(A.12) 

(A.13) 
(A. 14) 
(A.15) 



The mass-squared matrix of stops is given by 



Mi 



,2* 



m l2 m 22 



2 > 



where 



m 



11 



rn 



1(4 - al) (pI - pI) + \vIpI 



m 



22 



tL 8 V 3 

2 I 2/2 2 \ i 22 
m ^ + q9i{Pi-P2) + T ) VtP2i 



Vt 



12 



[(/ipi + AP2 cos 6 1 ) - iA t p 2 sin 5] . 



(A.16) 

(A.17) 

(A.18) 
(A.19) 



Here y t is the top Yukawa coupling. rn~^, m 2 R and A t are the soft-SUSY-breaking mass- 
squared parameters. 

The formulas for the corrections to the effective parameters are 



Apn 2 3 
A t ~A 6 



-N c y^A* t -i I A 1 (k)A 2 (k), 

Jk 



Jk 



(A.20) 
(A.21) 



-N c y^A*-i f 
Jk 



92 9t 



9\ 



2£-gJ A'i(k)A 2 (k) + ^A^AKk) 
+y^\ 2 Al(k)Al(k) 



(A.22) 



A,~A 7 = -N c y*fiA* t -i J . 



y 2_(9l_9T 
yt y 4 12 

2 1 A |2 A 2 



Al(k)A 2 (k) + ^-^\A 1 (k)Ai(k) 



+yi\A t \ 2 A((k)Ai(k) , 



(A.23) 



where 



Ai(A;) 



k 2 — mj ' 



A 2 {k) 



1 



k 2 



(A.24) 



tL tji 

0, the finite-temperature Feynman integrals become ill-defined 
because of infrared divergence. Then one should employ another method to evaluated the 



In the case of m tE ~ 



effective potential as done in [[[(J • B 2 and C\ are extracted from the light stop contribution 
to the effective potential, which is treated by the high-temperature expansion ( |2 . 1 3| ) . 
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